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Abstract 



Scalar perturbations can grow during a phantomic cosmological phase 
as the big rip is approached, in spite of the high accelerated expansion 
regime, if the equation of state is such that ^ = a < — |. It is shown 
that such result is independent of the spatial curvature. The perturbed 
equations are exactly solved for any value of the curvature parameter 
k and of the equation of state parameter a. Growing modes are found 
asymptotically under the condition a < — |. Since the Hubble radius 
^ ' decreases in a phantom universe, such result indicates that a phantom 

scenario may not survive longtime due to gravitational instability. 

(S| . PACS: 98.8G.-k, 95.36. +x 

' In the ordinary theory of cosmological perturbation [Ij El [3] there are two 

fundamental regimes defined from the notion of Jean's length, Xj. Perturbations 
whose scales are such that A > Aj suffer gravitational instability since the 
gravitational attraction dominates over the pressure reaction to contraction; on 
the other hand, perturbations whose scales satisfy the condition A < A j , tends to 
oscillates due to the effectiveness of the pressure opposition to the gravitational 
5J] ' collapse. In the relativistic version of the theory of cosmological perturbation we 

5^ I must add a new relevant scale, given by the Hubble radius. The Hubble radius 

define, in some sense, the effective causal region, and consequently the region 
where the effects of the microphysics phenomenon (responsible for the pressure) 
may play a relevant role in the process of gravitational collapse. Perturbations 
whose scales are much large than the Hubble radius tends to become frozen, 
while for those smaller than the Hubble radius, the pressure tends to produce 
a damping effect. Moreover, if the spatial section of the metric representing 
the universe is not flat a new scale appears which is connected to the curvature 
parameter. 



^e-mail: fabris@pq.cnpq.br. On leave of absence of Departamento de Fisica, Universidade 
Federal do Espi'rito Santo, 29060-900, Vitoria, Espirito Santo, Brazil 
■^e-mail: dfaragoObol. com.br 
•^e-mail: sergio.vitorino@pq.cnpq.br 



1 



In cosmology, concomitant with the interplay between the gravitational at- 
traction and the pressure resistance, there is a supplementary relevant effect 
due to the expansion of the universe, which acts as a friction term in the funda- 
mental equations of the perturbed quantities: the expansion leads to a damping 
in the evolution of perturbations. In general, for very large perturbations, that 
are connected with the large structures existing in the universe, there are two 
modes: a growing or a constant mode and a decreasing mode. However, if the 
pressure is negative, and the expansion becomes accelerated, gravity may be- 
come repulsive and we generally find only decreasing modes: the perturbations 
are always damped and structures can not be formed. 

There are quite strong evidences that the universe today is in an acceler- 
ated expansion phase [4]. If this is true, the material content of the universe 
must be dominated by an exotic fluid whose pressure is negative. There are 
many candidates for this exotic fluid, like cosmological constant, quintessence, 
K-essence, Chaplygin gas, etc., each of them having its advantages and disad- 
vantages from the theoretical point of view. We address the reader to recent 
reviews on the many dark energy models existing in the literature [51 [5] . Such 
negative pressure fluid may generate repulsive effects driving the accelerated 
expansion. There are some claims that the observational data favors a phan- 
tom fluid, that is, a fluid whose negative pressure violates the dominant energy 
condition p + p > 7J. Representing such fluid by a self- interacting scalar 
field, a phantom fluid requires a "wrong" sign in the kinetic term. This may 
imply instability at quantum level. However, interesting propositions have been 
made in what concerns such quite exotic fluid, mainly in the context of ghost 
condensation existing in some string configurations [9]. 

Classically, one of the most striking feature of phantom fluids is the fact that 
its energy density grows with the expansion. This is consequence of the violation 
of the strong energy condition, and it may lead to a future singularity in a finite 
future time, the so-called big rip. This is, of course, a very undesirable feature. 
But it has been shown [S] that in a single fiuid approximation, for a spatially flat 
universe, the scalar perturbations can grow when the scales of the perturbation 
are greater than the Hubble radius. Since, the Hubble radius decreases, for 
a phantom dominated universe, the isotropy and homogenous condition would 
not be satisfied anymore as the big rip is approached, leading perhaps to the 
avoidance of this future singularity, leaving a very inhomogeneous universe. This 
situation may occur if the pressure is negative enough in order to satisfy the 
condition - = a < — |. 

p o 

In the present work we extend those result showing that phenomenon of 
enhancing of the inhomogeneities in large scales is independent of the spatial 
curvature, and that the critical point a = — | is present in any class of homo- 
geneous and isotropic universe. In order to do so, we will solve the perturbed 
equations for scalar modes for any value of k and a. An asymptotic analysis 
will reveal the existence of critical behavior associated to a = — |. 

For a universe dominated by a fluid whose equation of state is given by 
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p = ap, the relevant equation is 
a'2 . SttG 



-pa 



p = Poa 



-3(l+a) 



(1) 



The primes indicate derivations with respect to the conformal time 77, and k is 
the spatial curvature parameter, k = ±1,0. The solution for this equation may 
be written in a unified form as 



a{v) = ao 



1 



Vk' 



1 + 3a 



Vkr 



(2) 



For k = 1, this solution represents a universe that begins and end at a singularity 
at a = for 1 < a < — ^, while for a < — ^ it represents a bouncing universe. 
On the other hand, for k — —1, 0, it represents an ever expanding universe for 
any value of a, with the following characteristics: for 1 > a > — ^, the expansion 
implies < 77 < 00 corresponding, in terms of the cosmic time i, to < t < 00; 
for — i > a, the expansion implies —00 < rj < 0, corresponding to < i < 00 if 



> a > — 1 and 



<i<Oif— l>a. This last feature leads to the notion 



of big rip. 

To study the evolution of scalar perturbations, we use the gauge invariant 
formalism. For a perfect fluid the evolution of the scalar perturbation is given 
by a single equation for the gravitational potential $ [3l [10] : 



$" + 3(l + a)i7$'- 



2H' + (1 + 3a){H^ -k)}<i> = 



(3) 



where H = ^ and 71^ is the eigenvalue of the Laplacian operator = —n?^. 
The perturbed equation can be recast under the following form 



(1 - z)z<S>" + 7^^i^(l - 2z)$' + = 



where 

1 + cos(Vfc6') 



2(1 + 3a) 



4fc 
(l + 3a)2 



(4) 



an^ -2{l + 3a)k 



1 + 3q: 



-V 



(5) 

The solution of Q can be represented under the form of hypergeometric 
functions. It reads in general, for any value of k and a, as follows: 



A+,A_-B-z 



where 



1 r„ 1 + a 



A+ - B + 1, - B + 1;2 - B:z 



l + 3a 



± J36 



(l + a)2 
(l + 3a)2 



4n2 



B = 



7 + 9a 
2(1 + 3a) 



(6) 



(7) 
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and c, c are constants. 

The asymptotic analysis for the flat case fc = was made in reference [8]. 
In this case, the hypergeometric's function reduces to Bessel's functions. It has 
been shown that one of the two modes remains constant, in the long wavelength 
limit, for any value of a. In the same limit, however, the other mode decreases 
when a > — I , but grows with time when a < — | . For a — both modes are 
constants. In the case 7^ 0, such analysis, in terms of large or small wavelength 
limit, is more involved since, contrarily to the flat case, we have a scale given by 
the Hubble radius and another scale given by the curvature. It becomes easier, 
in this sense, and for our purpose more relevant, to consider the behavior in the 
extremes of the time interval. In order to perform this analysis, we must take 
into account some convenient transformation properties of the hypergeometric 
functions. In special, the following transformations will be useful (see reference 

my- 

2F, {A, B- C- z) = r(C-^r)r(C-f) ^^^ {AB-A + B~C + l-\-z) + 
(1 - -A^C-B-C-A-B^l-l-z) ;(8) 

nC)T{B-A)^ . , . ^. . , . 0.1 



2F1 (A, B- C- z) = r(i3)r(C - A) (-1)^' \A,A + l~C-A+l-B, _^ 
nC)V{A-B), . ^^00,. ^.o , . ..1 



nA)nC - B) (-1)^^-^^^^ \B,B^l-C-,B^l-A,-] .(9) 
Hence, we have the following asymptotic behaviors: 



^ 2FM,B;C-z) ^ zi-^ , (10) 

T{C)T{C ~ A~ B) 
T{C - A)V{C - B)^ 
T{C)T{A + B-C) 



T{C)T{C ~ A~ B) , ,r A R 



T{A)T{B) ' ^^^^ 

+ r{A)r{c~By > ■ ^ ' 

Using these expressions, we can determine the behavior of the perturbations 
in the two different extremities of the time interval, for each value of k. 

• k = 1. In this case the conformal time interval is < 77 < tt?— for 

— ' — H-3a 

a > — i and < 1] < for a < — ^ . Using the asymptotic expressions 

written above, we find the following behaviors: for a > — ^, there is 
initially two decreasing modes and, as the universe approaches the big 
crunch a,t rj — j^^, there are a constant mode and a growing mode; for 
— i > q; > — I there is initially, during the contraction phase, a growing 
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mode and a constant mode, and as the scale factor diverges in the other 
asymptotic, there is a constant mode and a decreasing mode; for a = — |, 
both modes are constant at the beginning of the contraction phase and at 
the end of the expansion phase; for a < — | , there is a constant mode and 
a decreasing mode at the universe begins to contract, and there are two 
increasing modes as the universe approaches the big rip. 

• k = —1. The range of the conformal time is < 77 < cx) for a < — ^ and 
~oo < 77 < for a < — i. The open case is more involved because, in 
opposition to the closed universe, the asymptotic behavior of the modes 
depends on the scale of the perturbation. Let us consider the situation 
where the eigenvalues of the Laplacian operator is null. Hence, when 
a > — i there is initially and in future infinity two decreasing modes. 
For — 2 < a < — I there is initially two decreasing modes, but in the 
future infinity there is a constant mode besides a decreasing mode; the 
case a = — I differs from the preceding one by the fact that in future 
infinity both modes are constant. Finally, when a < — |, both modes 
are initially decreasing but they become growing modes as the big rip is 
approached. 

When a = in all cases, the same features observed in the flat universe are 
reproduced here, since for this particular equation of state the matter density 
scales as the curvature parameter in the Friedmann's equation. 

The main conclusion of the previous analysis is that a phantom cosmologi- 
cal scenario is highly unstable against scalar perturbations under the condition 
of an isotropic and homogeneous background universe if the pressure is nega- 
tive enough, that is - < — |: the scalar perturbations grows as the big rip is 
approached. As in tfie flat case, the Hubble radius shrinks with time in the 
phantomic case. This means that the large scale approximation becomes es- 
sentially valid asymptotically for all perturbation scales in the phantomic case. 
The analysis was made using a perfect fluid material content. However, at large 
scales, it is expected that a more fundamental representation, using for example, 
scalar fields, must give the same results as the perfect fluid case jl2| . 

A curious feature of the above results concerns the critical point a — — |. 
As far as we know, it does not correspond to any energy condition (contrarily to 
a = — i and a = —1). Also, it seems not related neither to the Hubble param- 
eter, nor to deceleration parameter, or even to the statefinder parameters jl3]. 
However, it seems to be a general critical point for the analysis of perturbation, 
that is not reflected in the kinematical quantities like those we have quoted 
above. A dynamical system analysis of the background does not reveal any 
particular feature at a = — | [2]: it is not a critical point for the background. 
Moreover, in the perturbed equation ^ there is no explicit special structure 
for that particular value of the parameter a. The nature of this critical point 
must still be cleared up. 
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